Influence of non-local exchange on RKKY interactions 
in III-V diluted magnetic semiconductors 
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The RKKY interaction between substitutional Mn local moments in GaAs is both spin-direction-dependent 
and spatially anisotropic. In this Letter we address the strength of these anisotropies using a semi- 
phenomenological tight-binding model which treats the hybridization between Mn d-orbitals and As p-orbitals 
perturbatively and accounts realistically for the non-local exchange interaction between their spins. We show 
that exchange non-locality, valence-band spin-orbit coupling, and band- structure anisotropy all play a role in de- 
termining the strength of both effects. We use these results to estimate the degree of ground-state magnetization 
suppression due to frustrating interactions between randomly located Mn ions. 

PACS numbers: 75.50.Pp, 75.30.Et, 75.30.Gw, 75.20.Hr 



The current interest in diluted magnetic semiconductors 
(DMS) is fueled by possible applications in spintronics and by 
basic-science issues generated by the interplay between disor- 
der, spin-orbit coupling, and magnetic order. We concentrate 
on the prototypical III-V DMS Gai-^Mn^As, which, once 
interstitial Mn ions have been eliminated, exhibits robust ho- 
mogeneous ferromagnetism 1 1 ] with critical temperatures T c 
above 160 K for x > 0.05. It is generally agreed that the 
substitutional Mn ions are in Mn 2+ valence states that have 
S = 5/2, L = local moments, and that exchange inter- 
actions with As neighbors allow the Mn moments to interact 
via valence-band holes |2]. The effective exchange interaction 
between Mn moments is spatially anisotropic and, because of 
spin-orbit interactions, also anisotropic in spin space. This 
Letter is motivated primarily by theoretical interest O, |4J] in 
the role of anisotropies in determining the character of the 
magnetic ground state but has implications for other aspects 
of (III,Mn)V DMS ferromagnetism. 

The theory of (III,Mn)V ferromagnetism has been devel- 
oped in several directions. A simple phenomenological ap- 
proach approximates the valence-band 
holes by a host-semiconductor Kohn-Luttinger envelope- 
function Hamiltonian and couples them to randomly located 
Mn spins by a local, isotropic exchange interaction J p( j. This 
leads to a semi-quantitative description of many transport and 
magnetic properties, particularly in the high-carrier-density, 
high-T c systems that are free of compensating Mn intersti- 
tials. However, it has led to conflicting conclusions on the 
importance of exchange anisotropy. The RKKY interaction 
obtained by Zarand and Janko 1 3] is highly anisotropic in spin 
space, i.e., it depends strongly on the orientation of two spins 
relative to their connecting vector, but it is spatially isotropic 
because it starts from a local hole-impurity exchange interac- 
tion and uses a spherical approximation for the bands. Us- 
ing a more realistic 6-band envelope-function Hamiltonian, 
Brey and Gomez-Santos 01 find that both spin and real space 
anisotropies are weak. Their conclusion, however, depends in 
part on their momentum-space cut-off 1 1 1 ] for the exchange 
interaction J p( j, i.e., on atomic-length-scale physics not de- 



scribed realistically in the envelope-function approach. First- 
principles calculations 1H2I1 do not have these limitations, but 
are hampered by their extreme sensitivity to the placement of 
unoccupied and occupied d-orbital energies relative to the va- 
lence and conduction bands. In this Letter we address ex- 
change anisotropy using a realistic tight-binding model that 
combines virtues of these two different approaches and es- 
timate the bulk magnetization suppression due to frustrating 
interactions between impurity moments. Based on our results 
we also suggest a possible route toward higher transition tem- 
peratures in (III,Mn)V ferromagnets. 



Our theory is based on a Slater-Koster 111 311 tight-binding 
model, and on a perturbative treatment of pd hybridization, in 
which the band electrons are integrated out to yield a spin- 
only model l3LI4.l5l ll4lll5lll6ll . A similar model has recently 
been used to obtain the local density of states around Mn im- 
purities 11711 . In Slater-Koster theory, the electronic structure 
is specified by orbital-dependent onsite energies and hopping 
amplitudes that are treated as fitting parameters. Spin-orbit 
coupling is included [18] to obtain realistic bands and a real- 
istic description of (III,Mn)V ferromagnetism |8]. 

Our Hamiltonian reads H = H c + Hd + i?hyb, where 



H B = 



k ( 



(!) 



describes perfect GaAs lll3ifl8il . Here, creates an elec- 
tron with wave vector k in orbital a with spin a. The most 
important effect of Mn impurities is to introduce partially 
filled d-orbitals. The resulting strong electron-electron in- 
teractions are parametrized by the local Hubbard repulsion 
U and the Hund's-first-rule coupling Jh fl9l EoTl : H c i = 
(e d + J H - U/2)N + 1/2 (U - J H /2)N 2 - J H S • S, with 

N = T,na d Ldna and S = f^nw* d na (<W / 2 ) d nv> , 

where d] ll7 creates an electron in d-orbital n with spin a. We 
assume U « 3.5 eV |H and J H ~ 0.55 eV |'22|. H hyh de- 
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scribes the hybridization between the d-orbitals and sp-bands, we can obtain analytic expressions from tight-binding theory. 
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where J\f is the number of unit cells in the system. The co- 
efficients are expressed in terms of real-space hopping matrix 
elements, t-^ an = J2i Ui Ucm, where the sum runs over 
nearest-neighbor As sites of the impurity. The symmetries of 
t\tan are obtained from Slater-Koster theory 1 13], which ex- 
presses the matrix elements in terms of two-center integrals. 
We use {pda) = 1.0 cV and (pdir) = —0.46 eV as inferred 
from photoemission l2lll and (sda) = 1.5 eV obtained as a 
rough spin average of ab-initio calculations for zinc -blende 
MnAs Si. 

In the large- U limit we can use canonical perturbation the- 
ory (CPT) [24] to integrate out d-shell charge fluctuations, 
leaving only the impurity spin degrees of freedom. We first 
consider a single Mn impurity. We introduce the canonically 
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transformed Hamiltonian H 
where T is hermitian, and expand in e. The operator T is cho- 
sen so that the linear term vanishes. To obtain manageable 
expressions we neglect the energetic spread of virtual band- 
electron states compared to the energy difference ~ JJ be- 
tween different Mn valence states. To be consistent we ignore 
contributions from bands other than the heavy-hole, light- 
hole, and split-off bands. Truncating the expansion at second 
order and projecting onto the N = 5, 5 = 5/2 ground-state 
subspace, we obtain 



H = H r 
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We have used that H^ h applied to a state in the (N, S) = 
(5, 5/2) sector results in a state with sharp quantum numbers 
(N, S) — (6, 2) and (4, 2), respectively. Ens is the corre- 
sponding isolated-ion energy. Inserting Eq. (|2j and noting that 
Zw d na (°w/ 2 ) d na' = S/5 in the (5, 5/2) sector, we ob- 
tain a Hamiltonian that includes a microscopic hole-impurity 
exchange interaction, 



H = H r 



(charge scattering) 
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The two energy denominators in 1/A are respectively the 
isolated-ion d 5 — > d A and d 5 — ► d 6 transition energies measured 
from the chemical potential. If either of the denominators be- 
comes small, the interval of energy over which our approx- 
imations are justified is correspondingly reduced. Note first 
that the exchange interaction is quite generally invariant under 
spin rotation. The wavevector dependence of the exchange in- 
teraction is specified by the factor ^ n tw a ' n for which 



For k, k' — > and a = a' = p x ,p y ,p z we obtain 
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3(pda) 2 - 4V3 {pda)(pdTr) + A{pdvf 



(6) 

Restoring the prefactor from Eq. (0J we find a microscopic 
expression for the envelope-function exchange constant J p d- 
By including the full (k, k') dependence we recover spatial 
anisotropies neglected in that theory. 

Since both denominators in 1/A must be negative for 
(5,5/2) to be the isolated-ion ground state, the exchange 
interaction is antiferromagnetic, J pc \ < 0. |J P d| is mini- 
mized and the effective model has the widest range of va- 
lidity when the d 5 — ► d 4 and d 5 — > d 6 transition energies 
bracket the Fermi energy Ep symmetrically. In this case 
J pc j = —48 meVnm 3 , close to the experimental value in 
(Ga,Mn)As 1 25 ] . We consider this case in what follows. The 
expression for J pc \, combined with materials trends l2fill . sug- 
gests that T c of Gai_ x Mn x Asi_ a Pj, quaternary alloys might 
increase with y since their d 5 — ► rf 4 transition energy will ap- 
proach Ep, increasing the value of J pc \. 

We employ the full (k, k')-dependent hole-impurity ex- 
change to evaluate the RKKY interaction between two Mn 
spins at and R and perform the CPT as above. Integrat- 
ing out the band electrons and expanding the action to second 
order in impurity spins we obtain 
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x {-no + c(k) - m)^ i M (k, k') (-iw + e(k') - ^y 1 

xj"(k',k) = -^V(R)^, (7) 

in/ 

where e(k) is the tight-binding Hamiltonian with ma- 
trix elements e a ' cr ' ;Q(T (k) and j' M (k, k') a v';aa = 
E n *Ln^a'n^ r The trace in Eq. Q is over orbital 
and spin indices. We diagonalize e(k) = d(k) U^, 
where d(k) is the diagonal matrix of band energies d Q(T (k), 
and perform the Matsubara sum. It is useful to express 
J M „(R) = / d 3 g/(27r) 3 e lq R J M „(q) in terms of its Fourier 
transform. Making use of the symmetries of d and U we 
obtain 

■^fa) = 9^2 / (27r1 3 E ^ kQCT E ^ _ /k-q,a'<r') 
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d Q v(k- q) - d a<T (k 
x [J7 k _ q ?(k-q,k)^] QV , ;Q(T 



f/ k ^(k,k-q)C/T_ q ] Q(T;QV , 
(8) 



where v uc is the unit-cell volume and f^ aa is a Fermi fac- 
tor. In the following, we take the electrons to be at T = 0. 
We remark that Eq. l|8) is unreliable when q is comparable 
to Brillouin-zone dimensions because the band eigenenergies 
are then as far from the Fermi energy as the d-quasiparticle 
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levels. Correspondingly the results for J^(IV) are quantita- 
tively reliable only for R 3> a, where a is the dimension of 
the fee unit cell. 
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FIG. 1: (color online) Fourier-transformed RKKY interaction 
Jxx(<l) and numerical errors in the (100) and (111) directions for 
Ef = —0.307 eV relative to the valence-band top, correponding to 
a hole concentration of 3.8 x 10 20 cm -3 . The dashed curve shows 
J xx (q) in the (100) direction calculated with a band Zeeman splitting 
of 0.125 eV, corresponding to 5% Mn substitution and full polariza- 
tion of Mn moments. 



We have evaluated J M „(q) using Monte Carlo (MC) in- 
tegration with the Vegas algorithm lE7il . Figure ^ shows 
Jxxil) in the (100) and (111) directions. At a nonzero Mn 
density, the interactions between spins are dominated by the 
pairwise RKKY interaction only if the mean hole-impurity ex- 
change interaction is weak Xji\$. This is indeed the case since 
Fig-fflshows that the effect of a realistic Zeeman splitting on 
J(q) is small. We note that J M „(q = 0) is isotropic; this 
limit determines the bulk magnetic anisotropy |9, 10] which 
vanishes in the present approximation lEill . 

Jpv (R) is evaluated as a Fourier sum over J M „ (q) calcu- 
lated on a cubic grid with (2nk) 3 /2 points in the fee Brillouin 
zone, making use of all symmetries. The resulting RKKY in- 
teraction is plotted in Fig. |2| It is ferromagnetic at small sep- 
arations, as expected. The near-neighbor interactions are not 
reliable, both because their evaluation stretches the validity of 
the CPT and because we neglect the superexchange interac- 
tion, which appears at fourth order in i?hyb, and in which an 
electron hops virtually from a Mn d-orbital to a d-orbital on 
a neighboring Mn site via an intervening As p-orbital. For 
larger separations J^ U (R) shows typical Friedel oscillations. 

We find a very strong anisotropy in real space; J p „(R) de- 
pends on the direction of R for similar R = |R|. This is 
a consequence of both the directionality associated with pd 
hybridization and of the anisotropy of the band structure and 
the Fermi surface; neither effect is included in the spherical 
model of Ref. |3J. In Ref. |4] the real-space anisotropy is con- 
cluded to be small, based on the interaction between two spins 
at neighboring sites. For small R we also find relatively weak 
anisotropies but at larger R this conclusion does not hold. The 
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FIG. 2: (color online) (a) Diagonal and (b) off-diagonal components 
of the RKKY interaction J^y (R) in various crystal directions, scaled 
by (R/a) 2 . All results have been obtained with rife = 36 and 2 x 10 s 
MC points for each q point except for (qa/2) 2 < 0.5, when 2 x 10 6 
points have been used. The off-diagonal components vanish exactly 
along (100). 



isotropic Gaussian ansatz for the hole-impurity exchange in- 
teraction 0] contributes to this small anisotropy. 

The anisotropy in spin space, i.e., the deviation of J M „(R) 
from J(R) S^u, is also large. For small spin-orbit coupling, 
the differences between diagonal components are of second 
order in spin-orbit coupling, whereas the off-diagonal com- 
ponents are linear. Only for the smallest separations is the 
relative anisotropy below 10% as found in Ref. |4|]. At larger 
R the anisotropy becomes quite pronounced, as in Ref. |3]. 

When the anisotropies are neglected, the moments are fully 
aligned in the ground state. To determine whether or not 
the anisotropies substantially alter the character of the ground 
state, we start from a fully aligned (in the z direction) spin 
configuration and consider the mean effective fields acting on 
individual spins, H^(Rj) = SJ^j^i Jfiz(R-i — Rj)> where 
the sum is over Mn impurity sites. Assumming that the Mn 
ions are distributed completely at random 1 16, 29], the aver- 
age over all sites is H ^ = (xS/v uc ) J^ z (q = 0) cx b^ z . On 
average the effective fields align with the average moment, but 
spatial fluctuations reduce the overall degree of spin polariza- 
tion. The typical angle of the Mn tilt at a given site is propor- 
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tional to the ratio of the xy plane effective-field components to 
H r . We find 



H 2 

-==?- = (x- 1 - l)v uc 



d 3 q |J, z (q)| 2 
(2tt)3 Jf,(q = 0)' 



(9) 



Thus the anisotropics become more important for small 
Mn fractions x. For the parameters used above we get 
~H2/[H Z ] 2 = 3.1 x lCT 5 ^- 1 - 1). We conclude that 
the anisotropics do not cause a large moment suppression in 
(Ga,Mn)As even for x ~ 0.01, despite the large anisotropics. 
The effect is small because many moments contribute to the 
effective field due to the long-range interaction, averaging 
over the anisotropics. We neglect the indirect influence of 
charge scattering, as well as Coulomb interactions and local 
chemical shifts. These will reduce the RKKY interaction at 
large separations and further reduce the importance of frus- 
trating interactions flUl . 

To conclude, we have used a Slater-Koster tight-binding 
model of III-V DMS to calculate the full momentum depen- 
dence of the hole-impurity exchange interaction. We find that 
this interaction depends crucially on the position of the Mn 
d-levels relative to the valence band and suggest that quater- 
nary Gai-^Mn^Asi-yPy alloys might have higher transition 
temperatures than Gai_ x Mn x As. Starting from the hole- 
impurity interaction, we have calculated the hole-mediated 
RKKY interaction between impurity spins. This interaction 
is highly anisotropic in real and spin space. The anisotropy 
crucially depends on two factors partly ignored in previous 
works: the nonlocal form of the hole-impurity exchange in- 
teraction and the highly anisotropic band structure. However, 
despite the strong anisotropies the local-moment suppression 
is weak due to the averaging brought about by the long-range 
RKKY interaction. 

We gratefully aknowledge stimulating discussions with W. 
A. Atkinson, L. Brey, T. Dietl, G. Fiete, T. Jungwirth, R Kac- 
man, T. Schulthess, J. Sinova, G. Zarand, and A. Zunger. 
AHM was supported by the Welch Foundation, by the Na- 
tional Science Foundation under grant DMR-01 15947 and by 
the DARPA SpinS program. 



* Electronic address: timm@physik.fu-beiiin.de 
t Electronic address: macd@physics.utexas.edu 
[1] S. J. Potashnik et al, Appl. Phys. Lett. 79, 1495 (2001); K. W. 
Edmonds et al, Phys. Rev. Lett. 92, 037201 (2004); D. Chiba 
et al, Appl. Phys. Lett. 82, 3020 (2003) and references therein. 
[2] H. Kepa et al., Phys. Rev. Lett. 91, 087205 (2003); Y. H. Mat- 
suda et al, cond-mat/0404635 (unpublished). 



[3] G. Zarand and B. Janko, Phys. Rev. Lett. 89, 047201 (2002); G. 
A. Fiete, G. Zarand, and K. Damle, Phys. Rev. Lett. 91, 097202 
(2003). 

[4] L. Brey and G. Gomez-Santos, Phys. Rev. B 68, 1 15206 (2003). 

[5] T. Dietl, A. Haury, and Y. Merle d'Aubigne, Phys. Rev. B 55, 
R3347 (1997). 

[6] T. Jungwirth et al, Phys. Rev. B 59, 9818 (1999). 

[7] T. Dietl et al, Science 287, 1019 (2000). 

[8] J. Konig, J. Schliemann, T. Jungwirth, and A. H. MacDon- 
ald, Electronic Structure and Magnetism of Complex Materials, 
Springer Series in Material Sciences 54, edited by D. J. Singh 
and D. A. Papaconstantopoulos (Springer, Berlin, 2003), p. 163. 

[9] M. Abolfath et al, Phys. Rev. B 63, 054418 (2001). 
[10] T. Dietl, H. Ohno, and F. Matsukura, Phys. Rev. B 63, 195205 
(2001). 

[11] J. Schliemann, J. Konig, and A. H. MacDonald, Phys. Rev. B 
64, 165201 (2001). 

[12] See for example S. Sanvito, P. Ordejon, and N. A. Hill, Phys. 
Rev. B 63, 165206 (2001); M. van Schilfgaarde and O. N. 
Mryasov, Phys. Rev. B 63, 233205 (2001); P. Mahadevan and 
A. Zunger, Phys. Rev. B 69, 115211(2004); L. M. Sandratskii, 
P. Bruno, and J. Kudrnovsky, Phys. Rev. B (to be published), 
cond-mat/0404163; T. C. Schulthess et al. (unpublished) and 
references therein. 

[13] J. C. Slater and G. F. Roster, Phys. Rev. 94, 1498 (1954). 

[14] G. Bouzerar, J. Kudrnovsky, and P. Bruno, Phys. Rev. B 68, 
205311 (2003). 

[15] D. J. Priour, E. H. Hwang, and S. Das Sarma, Phys. Rev. Lett. 

92, 117201 (2004); S. Das Sarma, E. H. Hwang, and D. J. 

Priour, cond-mat/0403059 (unpublished). 
[16] C. Timm, J. Phys.: Condensed Matter 15, R1865 (2003). 
[17] J.-M. Tang and M. E. Flatte, Phys. Rev. Lett. 92, 047201 (2004). 
[18] D. J. Chadi, Phys. Rev. B 16, 790 (1977). 
[19] R. H. Parmenter, Phys. Rev. B 8, 1273 (1973). 
[20] P. Kacman, Semicond. Sci. Technol. 16, R25 (2001). 
[21] J. Okabayashi et al, Phys. Rev. 59, R2486 (1999). 
[22] L. Craco, M. S. Laad, and E. Miiller-Hartmann, Phys. Rev. B 

68, 233310 (2003). 
[23] S. Sanvito and N. A. Hill, Phys. Rev. B 62, 15553 (2000). 
[24] J. R. Schrieffer and P. A. Wolff, Phys. Rev. 149, 491 (1966); 

K. A. Chao, J. Spalek, and A. M. Oles, Phys. Rev. B 18, 3453 

(1978). 

[25] H. Ohno, Science 281, 951 (1998); H. Ohno, J. Magn. Magn. 

Mat. 200, 110(1999). 
[26] A. Zunger in Solid State Physics 39, 275 (1986) and references 

therein. 

[27] G. P. Lepage, J. Comp. Phys. 27, 192 (1978). 

[28] The bulk magnetic anisotropy arises from weak multi-spin ef- 
fects neglected here that give rise to anisotropies consistent with 
cubic symmetry and from strain effects that lift the symmetry. 
Our focus here, however, is on frustrating spin interactions that 
could reduce the magnitude of the bulk magnetization. 

[29] C. Timm, F. Schafer, and F. von Oppen, Phys. Rev. Lett. 89, 
137201 (2002). 



